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Rapid AGN accretion from counter-rotating discs 



(N 

o 

(N 

X) 

tin 
On 



X 

Oh! 

o 



> 

00 

o 
o 
p 

en 

o 

(N 



t 



Christopher J. Nixon^*, Andrew R. King^ and Daniel J. Price^ 

^ Theoretical Astrophysics Group, University of Leicester, Leicester LEI 7RH UK 

^ Monash Centre for Astrophysics (MoCA), School of Mathematical Sciences, Monash University, Vic. 3800, Australia 



24 February 2013 



ABSTRACT 

Accretion in the nuclei of active galaxies may occur chaotically. This can produce accretion 
discs which are counter-rotating or strongly misaligned with respect to the spin of the cen- 
tral supermassive black hole (SMBH), or the axis of a close SMBH binary. Accordingly we 
consider the cancellation of angular momentum in accretion discs with a significant change 
of plane (tilt) between inner and outer parts. We estimate analytically the maximum accretion 
rate through such discs and compare this with the results of Smoothed Particle Hydrodynamics 
(SPH) simulations. These suggest that accretion rates on to supermassive black holes may be 
larger by factors > 100 if the disc is internally tilted in this way rather than planar. This off'ers 
a natural way of driving the rapid growth of supermassive black holes, and the coalescence of 
SMBH binaries. 

Key words: accretion, accretion discs - black hole physics - galaxies: formation - galaxies: 
active - cosmology: theory 



1 INTRODUCTION 

Astronomers still do not know how the supermassive black holes 
(SMBH) in galaxy centres accrete gas and grow. The problem lies 
in the very small angular momentum this gas must have in or- 
der to accrete on a reasonable timescale. For example, the viscous 
time of a standard accretion disc of ~ 1 pc radius accreting on 
to a SMBH of lO^M© approaches a Hubble time. A natural rem- 
edy for this is to assume that on small scales near the hole, accre- 
tion occurs in episodes whose angular momenta are randomly ori- 
ented with respect to each other. This chaotic picture predicts that 
the spin of the hole remains low, allowing r apid mass growth pro- 
vide d that there is an adequate mass supply ("King & Pringl e 200^ : 
[King &Pringle 2007: King et al. 2008 ; Fanidakis et al. 201l|). A 
characteristic feature of this type o f accretion is that the disc flow 
on to the SMBH is often retrograde (JKing et al.. 20051) with respect 
to the hole spin, or to the rotation of an SMBH binary resulting 
from a galaxy merger. 

In a recent paper (INixon et al., 20 lib we showed that a retro- 
grade coplanar external gas disc can be markedly more efficient in 
shrinking an SMBH binary than a prograde one. A retrograde disc 
can cancel the binary orbital angular momentum directly, rather 
than being slowed by the resonances which always arise in a pro- 
grade circumbinary disc. It is natural to ask if this angular momen- 
tum cancellation can work if we replace the SMBH binary with a 
pre-existing gas disc surrounding a single SMBH. If so, this may 
be a mechanism promoting much faster delivery of gas on to the 
black hole. This situation may arise naturally from misaligned ac- 
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cretion events, for example where an accretion event forms a disc 
around the black hole, and a second event forms an outer disc 
which is misaligned with respect to the first one. Similar cases oc- 
cur as a misaligned accretion disc attempts to reach an axisym- 
metric (co - or counter-aligned) configuration arou nd a spinning 
black hole JKing et al. 2005l : lLodato & Pringle 200^ . and also dur- 
ing the closely related co- or counter-alignment of an a ccretion 
disc around a central binary (INixon, King & Pringle, 20 lib . 

In all these cases integrating the full dynamics of the accre- 
tion disc is complex and time-consuming because the component 
discs can warp before interacting. To study the mixing of distinct 
or opposing disc angular momenta we investigate a simple form 
of this kind of accretion. We assume an accretion disc with initially 
distinct planes in its inner and outer parts. We let these two parts in- 
teract viscously and determine how much the accretion is enhanced 
compared with a disc lying in a single plane. In Section |2] we give 
simple analytical estimates of this. In Section[3]we perform global 
3D hydrodynamical simulations to examine the disc's dynamical 
behaviour in the light of these estimates. We interpret our results in 
SectionH 



2 COUNTER-ROTATING DISCS 

We consider an inner and outer disc counter-rotating with respect 
to each other but not coplanar, that is, the angle 6 between the disc 
angular momentum vectors obeys 7t/2 < 6 <n. 

If the discs were in perfect contact the gas would share op- 
posed angular momenta from almost the same radii. This may be 
unphysical (or imply a very rapid evolution away from the initial 
state) since the discs may attempt to open a gap of some size, as we 
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discuss later. For generality we choose an arbitrary gap between the 
discs of size a. 

Now we assume that the gas efficiently shares angular mo- 
mentum across this gap and that an equal mass of gas from each 
disc takes part in the interaction. The resultant angular momentum 
of the interacting gas is 
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L(/?gap + a) + L(/?gap)j 
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(1) 



to first order in a/Rg^p, where Rg^p is the outer edge of the inner 
disc (so that Rg^p + a is the inner edge of the outer disc), G is the 
gravitational constant, M is the mass of the central object, and we 
have assumed a Keplerian potential for the gas. 

The circularisation radius for this interacting gas is 
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We expect the size of the gap between the two discs to be roughly 
the scale of the epicyclic motions in the gas disc as this is where 
fluid orbits cross and interact. This implies a gap size ~ //(/?gap), 
where H is the semi-thickness of the disc, so that 
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For typical AGN disc parameters, H/R ~ 10~^, and there are two 
distinct cases. For 6 close to n, we can have 1 + cos 6 <^ H/R, so 
that 
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This suggests that the gas orbits would decrease by ~ 7 orders of 
magnitude in radius which implies direct accretion of the gas as 
long as its path to the hole is clear. 

If however 6 is not close to n we have 1 + cos 6 » H/R, and 
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Then the gas orbit decreases least when 6 ~ 7t/2, where Rdrc ^ 
\Rgap- The orbit scale 0]) varies smoothly between these two ex 
treme values. 

The accretion rate through a disc is approximately 

Mdisc 



M 

'^visc 

where Mdisc is the disc mass and TvIsc is the viscous timescale 
R' 



(6) 



(7) 



where Ris a. characteristic radius for the disc and v is a measure of 
the disc viscosity. Thus for constant viscosit3Q the accretion rate is 
oc R~^. This suggests that for 6 > n/1, where the gas orbits decrease 
by a factor > 4, there is potential for the accretion rate to increase 
by a factor > 16. For large ^ ^ :7r the cancellation can lead to direct 
accretion of gas on a dynamical timescale. This allows accretion on 
a local disc filling timescale at the radius of the gap. This timescale 
is much shorter than the viscous timescale as we only need move 
gas to /?gap and not right on to the hole. Note that the rate at which 
mass can be supplied to Rg^^ is clearly an upper bound on the pos- 
sible accretion rate from this process. 

^ This is only approximate as one expects y to increase with radius. 



If the gas does not accrete directly then the viscous timescale 
is shortened to 



l''^circ/-'^gap ) '^vi 



(8) 



Although these estimates are suggestive, we caution that we 
have assumed that the gas shares angular momentum efficiently, 
and that an equal mass of gas interacts from each disc. In reality 
neither of these simplifications may be valid. Accordingly we use 
Smoothed Particle Hydrodynamics (SPH) to simulate the interac- 
tions between the discs. 



3 SIMULATIONS 

3.1 Code Setup 

We use PHANTOM, a low-memory, highly efficient SPH code 
optimised for the study of non- self- gravitating problems. This 
code has performed well in related simulations. For example 
Lodato & Price (2010) performed simulations of warped accretion 
discs and foun d excellent agreement with the analytical work of 
lOgilvie (19991) on the nature of the internal accretion disc torques 
in response to warping. 

The implementation of accretion disc of-viscosity 
( Shakura & Sunvaev. 1973h in phantom is described in 
iLodato & Price (201Qh . Specificall y, we use the 'artificial v iscosity 
for a disc' described in Sec. 3.2.3 of lLodato & Price (2010h, si milar 
to earlier SPH accretion disc calculations (e.g. lMurrav. 19961) . The 
main diff'erences compared to standard SPH artificial viscosity are 
that the disc viscosity is applied to both approaching and receding 
particles and that no switches are used. Ou r implementation also 
diff'ers slightly from iLodato & Price (2010l) in that we retain the 
^^ term in the signal velocity in order to prevent particle inter- 
penetration from occurring in high- velocity stream-disc collisions. 
The disc viscosity in phanto m was extensively calib rated against a 
ID thin Qf-disc evolution in lLodato & Price (201Qh (c.f. Fig. 4 in 
that paper) and the disc scale heights employed here are similar. 
As the exact value of a is unimportant given the dynamical nature 
of the simulations, we simply use a^ - 1 and ^^ - 2 which 
at the employed resolution (see below) corresponds to a physical 
viscosity ofss ~ 0.02-0.06. 

For these simulations we model the gravity of the accretor 
with a Newtonian point mass potential, adopting an accretion ra- 
dius of 0.1 in code units. We start with a flat disc of gas, composed 
of 10 million SPH particles, in hydrostatic equilibrium between 1.0 
to 2.0 in radius, and surface density distribution 2 oc /?~\ setup 
using the usual Monte-Carlo technique. We also choose a vertical 
density profile corresponding to H/R = 0.02 at /? = 1 and employ 
an isothermal equation of state. To produce the two distinct planes 
for the disc we rotate the outer half (in radius) of the disc by an 
angle Q. We then perform simulations with diff'erent values of Q to 
find how the accretion rate changes with inclination angle. 



3.2 Tilted disc evolution 

We perform seven simulations corresponding to inclination angles 
Q - 30/1° where n - 0, 1, ...,6. The ^ = 0° simulation gives a flat- 
disc accretion rate to which we compare the other simulations. We 
note that the 180° case is a rather unrealistic setup as it would re- 
quire perfectly anti-parallel accretion events. Even in the co-and 
counter- aligned discs predicted by the numerical simulations of 
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Figure 1. The accretion rates for the isothermal simulations with different 
inclination angles 6. Time is in units of the dynamical time at 7? = 1 . The 
discs with a tilt > 90° are much more efficient at driving accretion, as these 
involve a degree of counter-rotation and so direct cancellation of angular 
momentum. 



iLodato & Pringle (2006h the discs never achieve a configuration 
where 6 is precisely 180°. 

The simulations were run with an isothermal equation of state, 
i.e. the temperature of each particle was held constant for each sim- 
ulation (see S ection [3 . 2 . 2 1 belo w for a discussion of the effects of 
different equations of state). This implies that the pressure, deter- 
mined hy P - cIp, is simply proportional to p, and the sound speed 
Cs is constant in time and the same for all particles. 

In Fig.[T]we report the accretion rates achieved for the different 
values of 6. This shows that for discs inclined by less than 90° the 
accretion rates are similar to that of a flat disc for the duration of our 
simulations. In this case the interaction between inclined gas orbits 
is relatively weak with the gas staying on near circular orbits, with 
planes slowly varying from the inner to outer disc. In contrast the 
simulations with 6 > 90° show a dramatic increase in the accretion 
rate. This is caused by the direct cancellation of orbital angular 
momentum and energy where the orbits of the disc gas cross. 

In Fig. [2] we show the disc structure for the simulation where 
6 = 30°. Here the two discs initially join to form a coherent disc 
with a warped region joining the still misaligned inner and outer 
discs. Gas is quickly depleted from the warp region, leaving a clear 
drop (of 2-3 orders of magnitude) in projected density in the warp 
region. This break between the discs persists throughout the simula- 
tions as the two discs slowly align with each other, reaching an an- 
gular separation of ~ 10° by the time the simulation ends. The large 
gradient in specific angular momentum in the warp region causes 
gas to quickly move on past the warp. T his causes the low dens ity in 
the warp region which is the disc break (INixon & King, 2012h . The 
disc therefore maintains two distinct planes. We still resolve the 
rings of gas between the two discs but thi s is clearly of much lo wer 
density than the inner and outer discs fcf iLodato & Price 20101) . In 
Fig. [3] we show the similar disc structure for the simulation where 

e = 60°. 

The 6 = 90° simulation (not illustrated) has an initial period 
of chaotic flow between the two discs, until the inner disc moves 
towards an inclination of ~ 85° after mixing with some of the gas 
from the outer disc. At this point a strong warp is set up and the 
flow continues in a similar way to ^ = 60° . 




Figure 2. Column density projection showing the disc structure for the 
6 = 30° simulation. The disc has evolved for ~ 50 dynamical times (-10 
orbits). 




Figure 3. Column density projection showing the disc structure for the 
6 = 60° simulation. The disc has evolved for ~ 50 dynamical times (-10 
orbits). 



The evolution of the 6 = 120° and 6 = 150° simulations pro- 
ceed in a dramatically distinct way to those with smaller 6. This is 
simply because the sharing of angular momentum now causes sig- 
nificant changes to the gas orbits. In the 6 = 120° simulation gas 
falls from the region between the two discs and circularises at a ra- 
dius inside the inner accretion disc. In this process some of gas is 
also accreted, as it carries a spread of angular momentum. Clearly 
this gas does not have precisely the angular momentum of the inner 
disc; it is a mix of gas from both discs. Therefore the quasi-steady 
state for this simulation is an outer disc causing cancellation with 
the outer edge of the original inner disc. This gas falls and feeds 
both direct accretion and the new innermost disc (white region of 
high density in Fig.|4]). This innermost disc is strongly warped with 
respect to the original inner disc and so there is again extra dissi- 
pation between these two discs as they try to straighten. The disc 
structure can be seen in Fig. |4] 

In the 6 = 150° simulation the gas falls to a much smaller ra- 
dius from the region between the two discs. Initially a large amount 
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Figure 4. Column density projection showing the disc structure for the 6 = 
120° simulation. The disc has evolved for ~ 50 dynamical times (-10 
orbits). Note the new innermost disc (of high density and therefore white) 
which has formed from the dynamical infall of gas from the warp region. 




Figure 5. Column density projection showing the disc structure for the 6 = 
150° simulation. The disc has evolved for ~ 50 dynamical times (-10 
orbits). 



of gas is directly accreted on to the sink particle representing the 
central accretor. The strong dissipation occurring early in the sim- 
ulation also causes the inner disc to spread all the way in to the ac- 
cretion radius. Some of the falling gas tries to circularise at a radius 
where it is forced to impact upon the inner disc. As it still has the 
sense of angular momentum of the outer disc this causes more can- 
cellation of angular momentum, driving more accretion. For these 
reasons this simulation generates the highest accretion rates - this 
is simply because it causes the greatest mixing of angular momen- 
tum in the gas. The disc structure for this simulation can be seen in 

Fig. a 

3.2.1 6* =180° 

The 6 = 180° case does not initially appear a very realistic ini- 
tial setup. However something close to this configuration may ap- 
pear in the evolution of accretion discs around spinning black 



Figure 6. Column density projection showing the disc structure for the 6 = 
180° simulation. The disc has evolved for - 50 dynamical times (-10 
orbits). The innermost disc is formed of gas falling from the gap between 
the original inner disc (now the middle disc) and the outer disc. In this 
picture the innermost and outermost discs are rotating clockwise and the 
middle disc rotates counter-clockwise. 



holes fcf iKing et al. 20051 : lEodato & Pringle 200d) and also dur- 
ing the evolution of an accretion disc around a binary system (cf 
iNixonetal. 201lh . 

Initially while the two discs (inner and outer) are in contact 
there is great cancellation of angular momentum for the gas. This 
gas then attempts to fall towards the central sink particle. It there- 
fore impacts upon the inner disc. The pressure of the inner disc gas 
is enough to force the falling gas out of the plane of the disc where 
it can continue on a low angular momentum orbit towards the cen- 
tre. It then circularises inside the inner disc (with some accretion) 
after colliding with gas on similar orbits and losing orbital energy 
in shocks. This new innermost disc clearly has the same sense of 
angular momentum as the outer disc since this has the larger spe- 
cific angular momentum. Once the original inner disc has spread 
inwards (and this new innermost disc has spread outwards), exactly 
the same process starts again between these two discs. See Fig. [6l 
for the disc structure. 

Angular momentum cancellation thus proceeds as follows. 
Imagine the two discs are not in contact. They spread viscously 
until gas from the discs can interact, sharing and so destroying an- 
gular momentum. Now consider two test particles (one from the 
inner disc and one from the outer disc) which interact. The par- 
ticle from the inner disc adopts a slightly eccentric orbit, which 
passes through its parent disc. It rejoins its own disc, which shrinks 
slightly to accommodate the loss of angular momentum. The parti- 
cle from the outer disc also adopts a slightly eccentric orbit (inside 
its original orbit). But this forces it to interact more with the 'hos- 
tile' inner disc. This causes runaway angular momentum loss for 
the outer disc particle, almost to the point of free-fall to the centre. 
Now the particle's free-fall energy is enough to move it out of the 
disc plane (helped by the pressure force from the inner disc) . It thus 
moves over and then inside the inner disc. The overall eff'ect is to 
shrink the inner disc and move the outer disc particles on to orbits 
that pass inside the inner disc. 

At this point one of three things can happen to the outer disc 
particle: 
1) it may have cancelled enough angular momentum to accrete di- 
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Figure 7. Column density projection showing the disc structure for the 6 = 
150° simulation with an isothermal equation of state. The disc has evolved 
for ~ 20 dynamical times. 



rectly. 

2) it may have cancelled enough to fall inside the inner edge of the 
discs where it crosses the plane of the inner disc and meets other 
particles on similar orbits: the particles shock and dissipate energy, 
forming a new disc at much smaller radii. This disc should be pla- 
nar, but any small perturbations in the initial discs will be amplified 
here at smaller radii and could produce a disc with a substantial 
warp. 

3) The particle's new orbit may take it through the inner (initial) 
disc. If this is the case it again interacts with the (hostile) disc and 
falls further. There are again two possibilities here. Either the par- 
ticle bounces all the way in past the hostile disc, or it has so many 
interactions that it adopts the sense of rotation of the hostile disc. 
So the particle either accretes or drives the inner disc to shrink, and 
hence drives accretion through the inner disc. 

All sequences here greatly enhance the central accretion rate 
while the discs are trying to mix their angular momenta. 



3.2.2 Effect of the gas equation of state 

We have also investigated the effect of changing the thermal treat- 
ment of the gas. In the simulations detailed above we used an 
isothermal equation of state. The gas was assumed to radiate away 
any dissipative heating instantly. The other extreme is to assume 
that any dissipation feeds the internal energy of the particles, which 
are not allowed to cool by radiative losses. Therefore we also simu- 
late the 6 - 150° case using an adiabatic equation of state (with 
7 = 5/3), transferring the kinetic energy dissipated by the vis- 
cous terms into internal energy. In all other ways this simulation is 
precisely the same as the corresponding isothermal simulation. In 
Figs.[7]&[8]we show the disc evolution after ~ 20 dynamical times. 
In the isothermal case the gas cancels angular momentum and falls 
to circularise at a new radius. In comparison Fig. |8] shows that the 
equation of state can play a significant role in the hydrodynamics. 
In this case the gas is significantly heated when it shocks and thus 
initially expands violently. This process drives more mixing of the 
gas and thus the accretion is now almost a factor of two higher here 
than in the isothermal case. The gas is still bound to the black hole 
and so returns on chaotic orbits to eventually accrete. 




Figure 8. Column density projection showing the disc structure for the 6 = 
150° simulation with an "adiabatic" equation of state (see the text for a full 
definition). The disc has evolved for ~ 20 dynamical times. Note that the 
gas no longer can cool to form a smaller disc (cf fig|7], but instead forms 
outflows which return on chaotic orbits generating a lot of cancellation of 
angular momentum. 



4 DISCUSSION 

Retrograde circumb inary accretion discs a r e natural in the chaoti c 
accretion picture dKing & Pringle 20061 : iKing & Prinde 20071) . 
and offer the possibility of ca ncelling the binary orbital angular mo- 
mentum (INixon et al.. 201lh . We have performed numerical simu- 
lations of accretion discs with two initially distinct planes. In the 
region between the two discs angular momentum is shared, driving 
accretion, particularly if ^ > ;7r/2. 

The highest accretion rates come from discs with inclination 
6 ^ 150°. Here the gas cancels angular momentum when the discs 
meet, allowing this gas to fall directly to the centre without the 
need to interact with the inner disc. Clearly this gas would in re- 
ality have some residual angular momentum and so would circu- 
larise at a new smaller radius. In the simulations however the gas 
falls inside our accretion radius for the central sink particle and so 
is directly accreted. Gas falling from between the discs always has 
angular momentum with the same sense of direction as the outer 
disc. Therefore if the gas tries to circularise outside the accretion 
radius it impacts the inner disc, cancelling more angular momen- 
tum and falling further until it is accreted. 

The accretion rates found here are > 100 times those of a flat 
planar disc. This increase comes from the mixing of angular mo- 
menta and consequent dynamical infall of gas to smaller radii. As 
gas is depleted from the gap the disc tries to spread and fill the 
gap, promoting more cancellation of angular momentum and so 
more accretion. The timescale for gas accretion is shortened (cf. 
equation[8]). If enough angular momentum is cancelled then the gas 
accretes directly on a dynamical timescale. 

Our simulations show that discs with an initial inclination 
less than 90° evolve into coherent discs joined by a warp. Over 
time, dissipation from this warp brings the two discs into a com- 
mon plane. This extra dissipation enhances accretion, the effect 
increasing with inclination angle (although negligible in compar- 
ison to the enhanced accretion when 6 > n/l). The discs appear 
to break, in the sense that the surface density of the disc is greatly 
reduced in the warp. However in these high-resolution simulations 
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we still resolve the orbits of the disc particles in the warp. This 
disc breaking is predicted by the numerical simulation s of warp 
propagation in accretion discs by Nixon & King (20121) using the 
constrained viscosities of lOgilvie ri999h. and has been also been 
seen in the SPH simul ations of lLarwood & Papaloizou (19971) and 
iLodato & Price (201Qh . 

Inclinations greater than 90° cancel rather than sharing angu- 
lar momentum between particles. This leads to particles falling on 
to smaller orbits, and so the accretion rates in these simulations are 
significantly higher. The accretion rate for 6 = 180° is lower than 
that of ^ = 150°. This is probably because gas must be pushed out 
of the plane of the discs before it can leap over the inner disc. 

The numerical simulations do not predict accretion rates as 
high as the analytical approach (|3]). This was based on two assump- 
tions: that the angular momentum cancellation is perfect, and an 
equal mass from each disc interacts. The simulations suggest in- 
stead that the gas falling from the gap has a spread of angular mo- 
mentum, where some accretes and some falls to a smaller radius. 

The main result of this paper is that internally tilted accre- 
tion discs can generate enhanced accretion rates of up to ~ 10"^ 
times that of a plane disc (cf Fig.[T]). In a quasi-steady state (such 
discs never reach a full steady state) we see accretion rates that are 
> 100 times those of a planar disc. As remarked above, discs tilted 
in this way arise naturally during the align ment of a disc around a 
sp inning black hole or a central binary (cf. iLodato & Pringle 20061 
& lNixon et al. 201 ih . We also expect such structures to occur dur- 
ing chaotic gas infall into galactic centres. This is a very efficient 
way of driving gas right down to the vicinity of an SMBH or an 
SMBH binary and thus causing accretion. Nested discs suggest a 
way of transferring gas from large radii (approaching galaxy scales) 
to small radii. 

We note finally that our result, that tilted and nested discs 
imply enhanced accretion, has implications for simulations where 
subgrid models are used to mimic feedback from the SMBH. The 
dynamics of misaligned accretion events may imply much larger 
accretion rates than predicted by the simple assumption of a planar 
disc with a fixed direction of angular momentum. 
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